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Abstract 

We address the problem of deconvolving the effects of interplanetary transport on observed 
intensity and anisotropy profiles of solar energetic particles with the goal of determining the time 
profile and spectrum of particle injection near the Sun as well as the interplanetary scattering 
mean free path. Semi-automated techniques have been developed to quantitatively determine 
the best fit injection profile, assuming (1) a general piecewise linear profile or (2) a Reid profile 
of the form [C/(t — t Q )] exp[— A/(t — t ) — (t — to)/B}. The two assumptions for the form of the 
injection profile yielded similar results when we tested the techniques using ISEE 3 proton data 
from the solar flare events of July 20, 1981 (gradual flare), and January 2, 1982 (impulsive flare). 
For the former event, the duration of injection was much shorter for protons of higher energy 
(75-147 MeV), which may be interpreted as indicating that the coronal mass ejection-driven 
shock lost the ability to accelerate protons to ~100 MeV after traveling beyond a certain 
distance from the Sun. 
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1. Introduction 

One of the basic questions one can ask about en- 
ergetic particles produced by solar activity is how 
many are accelerated as a function of energy and 
time. Electromagnetic and neutron diagnostics of in- 
teracting energetic particles can provide a wealth of 
information on the particle spectrum and the timing 
of particle acceleration at the flare site [e.g., Chupp, 
1990; Aschwanden et al, 1994; Aschwanden, 1996; 
Debrunner et al., 1997; Rank et al, 1997a, b]. Com- 
plementary information can be obtained from space- 
craft or ground-based observations of escaping, ener- 
getic, charged particles. The interpretation of these 
direct observations is complicated by the effects of 
interplanetary transport, particularly scattering due 
to irregularities in the interplanetary magnetic field 
[Meyer et al, 1956]. To date, such observations have 
always been performed at distances approximately 0.3 
AU or greater from the Sun, while the mean-free path 
of interplanetary scattering parallel to the magnetic 
field, Am, is typically between 0.08 and 0.3 AU (near 
the Earth), though it can be > 1AU for infrequent 
"scatter-free" events [Palmer, 1982]. Since the dis- 
tance from the Sun to the observer is of the same order 
of magnitude as, and usually larger than, the mean- 
free path, the interplanetary transport is largely diffu- 
sive in most cases, making it difficult to determine the 
underlying time profile of injection near the Sun and 
possibly affecting the measurement of the spectrum. 

Why should we go through the trouble of decon- 
volving the effects of interplanetary transport to de- 
termine the injection profile of escaping energetic par- 
ticles, when electromagnetic or neutron diagnostics of 
interacting particles can yield such information more 
directly? A key motivation is that it has recently be- 
come widely accepted that for many flare events the 
interacting and escaping energetic particles are accel- 
erated in different locations and possibly by differ- 
ent mechanisms. Twenty years ago, it was recognized 
from X ray observations that solar flares can be ba- 
sically classified into two groups [Pallavicini et al, 
1977], which are now commonly called "impulsive" 
and "gradual" according to the duration of their X 
ray emission. More recently, flares in these two cate- 
gories have been found to exhibit many other differ- 
ences in their physical properties and particle emis- 
sions, although measurements of ionic charge states 
indicated that ions from both types of flares are ac- 
celerated out of coronal material [Luhn et al, 1984, 
1987]. One key difference is that gradual flares are 
much more likely to be associated with coronal mass 



ejections (CMEs), which in turn often drive a trav- 
eling interplanetary shock [Sheeley et al, 1983]. For 
over a decade there has been mounting circumstantial 
evidence that for flare/CME events, the flare is re- 
sponsible for accelerating interacting particles, while 
the CME and/or associated shock is responsible for 
accelerating escaping ions on open magnetic field lines 
over a wide range of hcliolongitudcs [e.g., Mason et 
al, 1984; Lee and Ryan, 1986; Reames, 1990]. (The 
origin of escaping electrons is subject to debate, but 
there is some evidence that they come from the flare 
site [Droge et al, 1990a].) Further evidence against 
the alternative and previously accepted view that es- 
caping ions are accelerated at the flare site followed 
by transport within the corona is provided by recent 
Solar, Anomalous, and Magnetospheric Particle Ex- 
plorer (SAMPEX) measurements of the charge states 
of escaping ions [Leske et al, 1995; Mason et al, 1995; 
Oetliker et al, 1997], which rule out transport within 
the corona for longer than ~1 min [Ruffolo, 1997]. 
Therefore observations of escaping ions can provide 
unique information about the acceleration of ions out 
of coronal material in the CME/shock region [Kahler 
et al, 1990]. 

For impulsive flare events, the remarkable enhance- 
ment of the 3 Hc/ 4 Hc ratio among escaping ions [Hsieh 
and Simpson, 1970] by up to a factor of thousands is 
best explained by turbulent acceleration in the flare 
region [Fisk, 1978]. The observed longitudinal dis- 
tribution of solar energetic particles from impulsive 
events (albeit narrower than for gradual events) could 
be explained by interplanetary transport perpendic- 
ular to the average magnetic field or by transport 
within the corona. Thus, for impulsive events, a de- 
convolution of the effects of interplanetary transport 
along the field could be useful for characterizing the 
process of azimuthal transport, especially if one stud- 
ies data from multiple spacecraft at different loca- 
tions. 

Previous analyses of the intensity and/or aniso- 
tropy of solar energetic particles versus energy and 
time for various solar flare/CME events have yielded 
information on (1) the injection of particles near the 
Sun versus energy and time and (2) the mean-free 
path of interplanetary scattering. In most cases, both 
must be determined to provide a good fit to the data, 
although some previous studies have reported ade- 
quate fits by assuming an instantaneous injection of 
particles for some events. 

Previous fits to solar energetic particle observa- 
tions have employed various assumptions about in- 
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terplanetary transport. Some authors have restricted 
their analyses to cases where the mean-free path is 
believed to be long and have neglected interplanetary 
transport effects altogether [e.g., Debrunner et al., 
1988; Kahler et al, 1990; Kallenrode and Wibberenz, 
1991]. Analytic solutions employing the diffusive ap- 
proximation [Parker, 1963] were used by Wibberenz et 
al. [1989], and approximate analytic solutions of an 
equation of focused transport [Earl, 1976] were used 
for a comprehensive survey of solar energetic parti- 
cle observations by Ma Sung and Earl [1978]. Lock- 
wood et al. [1982] used a Monte Carlo method to solve 
the hard-sphere pitch angle scattering model of Fisk 
and Axford [1969]. Numerical solutions of a diffu- 
sive model that incorporates solar wind convection 
and adiabatic deceleration [Hamilton, 1977] were also 
employed by Beeck and Wibberenz [1986] and Beeck 
et al. [1987]. Numerical solutions of the equation of 
focused transport of Earl [1976] by an cigenfunction 
expansion technique were used by Bieber et al. [1980, 
1986], and solutions using a finite difference method 
[Ng and Wong, 1979; Schliiter, 1985] have also been 
employed [e.g., Droge et al., 1990b; Kallenrode et at, 
1992]. 

Naturally, a deconvolution that employs a less ac- 
curate transport model should yield a less accurate 
injection function. In particular, our results indi- 
cate that deconvolution with a less accurate transport 
model yields an artificially broad or "defocused" in- 
jection function. This underlines the importance of 
accurate modeling of interplanetary transport effects. 
A recent comparison [Earl et at, 1995] showed that 
nearly identical results were obtained by three inde- 
pendent computer codes for treating interplanetary 
scattering, which used a Monte Carlo method [Earl, 
1987], a finite difference method [Ruffolo, 1991], and 
an eigenfunction expansion [Pauls and Burger, 1994]; 
the consistency gives one confidence in the numerical 
accuracy of such methods. The transport equation 
and finite difference code of Ruffolo [1995] also include 
convection and adiabatic deceleration in the frame- 
work of focused transport. Solutions of this trans- 
port equation by other finite difference codes [Hatzky, 
1996; Lario, 1997] and by a Monte Carlo code (L. 
Kocharov, private communication, 1997) have cor- 
roborated these results. Again, the agreement makes 
one confident that the various methods are accurately 
solving the transport equation. Remaining system- 
atic errors could arise from the underlying transport 
assumptions. 

The goal of the present work is to apply such state- 



of-the-art transport simulations and develop semi-au- 
tomated fitting techniques to accurately determine 
the injection of particles near the Sun as a function 
of energy and time. The fits are objective, relying on 
X 2 minimization instead of eyeball evaluation as in 
much previous work. Two deconvolution techniques 
have been tested for protons from the gradual flare 
event of July 20, 1981, and the impulsive flare event 
of January 2, 1982. We have successfully fit the data, 
and the two techniques yield consistent results. For 
the gradual event of July 20, 1981, the duration of in- 
jection was much shorter for protons of higher energy 
(^100 McV), which we interpret as indication that 
the CME shock no longer accelerated protons to such 
high energies after traveling beyond a certain distance 
from the Sun. 

2. Numerical Techniques 

To simulate the interplanetary transport of solar 
energetic particles, we solve a Fokker-Planck equa- 
tion of focused transport that includes the effects of 
interplanetary scattering, adiabatic focusing, and so- 
lar wind effects, such as solar wind convection and 
adiabatic deceleration, to first order in (v sw /c), where 
v sw is the (constant) solar wind speed [Ruffolo, 1995]: 
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where F(t, /i, z,p) = d 3 N / (dzdfidp) , the density of 
particles in a given magnetic flux tube (following Ng 
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and Wong [1979]) as a function of the four indepen- 
dent variables: t (time), /i (pitch angle cosine in the 
solar wind frame), z (distance along the magnetic 
field), and p (momentum in the solar wind frame). 
Also, v is the particle speed, ip{ z ) is the "garden hose 
angle" between the magnetic field and the radial di- 
rection, L(z) is the focusing length, B/(dB/dz), ip(fi) 
is the pitch angle diffusion coefficient, and E' / E = 
1 — nvswvseap/c 2 is the ratio of the total energy in 
the solar wind frame to that in the fixed frame. 

We assume an Archimedean spiral magnetic field 
[Parker, 1958] for the observed solar wind velocity. 
The pitch angle scattering coefficient, ip, is parame- 
terized as 

^)=^r i a-M 2 ), 

following Jokipii [1971]. We have used q — 1.5, which 
is in the range of 1.3-1.7 inferred by Bieber et al. 
[1986]. The amplitude, A was determined from Am, 
the mean- free path parallel to the field, or A r , the 
radial mean free path, which are related by 

A r 3 v 

cos 2 V = 11 = (2 -q)(4- q) ~A' 

The initial condition places all the particle density at 
the inner boundary at r = 0.01 AU, simulating an in- 
stantaneous injection near the Sun. Thus the simula- 
tion results form a Green's function for the response 
to a J-function injection. The initial spectrum was 
estimated from the data. We use absorbing (zero in- 
flow) boundary conditions at the inner boundary and 
at an outer boundary far enough to be inaccessible to 
particles during the course of the simulation. 

We numerically solved (1) using the computer code 
of Ruffolo [1995], as modified to use t instead of vt 
as an independent variable because data gaps can 
make it awkward to bin the data according to vt on 
a particle-by-particle basis. Simulations using t or vt 
as an independent variable yielded consistent results. 
Each simulation over the time of interest for this work 
required ~4 hours on a Sun Ultra- 1 workstation. 

Since F is defined with respect to fi and p in 
the local solar wind frame, it is necessary to trans- 
form F into the fixed frame to predict counting rates. 
This is known as the Compton-Getting transforma- 
tion [Compton and Getting, 1935]. We have fit data 
from the MEH instrument [Meyer and Evenson, 1978; 
Kroeger, 1986] on board the ISEE 3 spacecraft, which 
are collected in eight orientational sectors as the in- 
strument's field of view (half-opening angle of 25°) 
rotates in the ecliptic plane. (Because of the narrow 



field of view of the MEH instrument, the measure- 
ments are essentially restricted to the ecliptic plane.) 
Therefore we used a Monte Carlo simulation to per- 
form the Compton-Getting transformation and cal- 
culate a matrix for converting F(fi,p) into predicted 
counting rates in the eight sectors at the energy of 
interest, taking into account the geometry of the de- 
tector. The simulation also took into account the 
predominant magnetic field direction during the time 
of interest, rotating the angular distribution so that 
/i = 1 pointed along that direction. 

The simulated and observed count rates have been 
compared in terms of the intensity and the aniso- 
tropy times intensity. The intensity is simply the 
sum of the eight sectored rates. We use the aniso- 
tropy times intensity instead of the anisotropy alone 
because the product can be approximated by a linear 
combination of sectored rates, which is necessary for 
one of our deconvolution techniques. In general, the 
observed magnetic field direction and the axis of sym- 
metry of the particle distribution both fluctuate with 
time, and the two do not track each other [Bieber and 
Evenson, 1987]. Presumably the axis of symmetry of 
the particle distribution is instead tracking a spatial 
average of the magnetic field over the particle gyra- 
tions. This makes it difficult to precisely predict the 
direction of the anisotropy, and therefore we compare 
the predicted and observed anisotropy times intensity 
values that are calculated with respect to the sector 
with the highest counting rate. Setting 9 = along 
that sector, a first-order harmonic expansion gives 
F = (F)(l + ScosO), where 9 is the angle and S is 
the two-dimensional anisotropy in the ecliptic plane. 
Then it can readily be shown that 5 = 2 (cos 9), where 
the average is weighted by the particle distribution, 
and the anisotropy times intensity is approximated by 
2 V Fj cos 9j , where Fj is the counting rate in sector 
j. (If the measurements were evenly distributed in 
three dimensions instead of two, it would be approx- 
imated by 3 ^2 j Fj cos#j.) Note that the anisotropy 
times intensity is calculated in the same manner for 
both the simulation results and the observations, so 
an error in the approximation should not strongly af- 
fect the comparison. 

We have developed two techniques for deconvolv- 
ing the effects of interplanetary transport in order to 
determine the underlying time profile of injection near 
the Sun, which in turn yields the injected spectrum 
as well as the best fit value of the interplanetary scat- 
tering mean- free path. These techniques solve the 
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inversion problem 

R(t) = f G(t- t')I{t')dt' , (2) 
Jo 

where I(t') is the injection of particles versus time 
near the Sun, G(t — t') is the Green's function, or 
the response to a ^-function injection, which is cal- 
culated by the transport simulation, and R(t) is the 
"response," i.e., the measured intensity or anisotropy 
times intensity at the spacecraft. Both techniques are 
semi-automated in that a computer program finds in- 
jection parameter values that minimize the x 2 of the 
fit between R(t) and the observational data (minus 
background) . Results for different values of A are then 
compared, and the overall best fit is used as the final 
result. 

The first deconvolution technique finds the optimal 
piecewise linear injection function for a given set of 
"joints," t.- L . To illustrate the technique, we consider 
fitting the intensity of 27-147 MeV protons measured 
by 1SEE 3/MEH after the gradual flare event of July 
20, 1981 (Figure la). The times of the joints in the 
piecewise linear injection function, I(t'), are chosen 
a priori, and the injection function is constrained to 
be zero at the first and final joints (Figure 2c). Then 
I(t') is a linear combination of triangular injections, 
Ii(t') (Figure 2a). The first triangular function starts 
from at the first joint, rises linearly to 1 at the next 
joint, and falls linearly to at the following joint (in 
units of 10 26 si-" 1 s" 1 MeV -1 , i.e., 10 26 per unit solid 
angle of the solar surface, time, and energy). The 
peak time of this function is then the start time of the 
next function and so on. We then convolve G(t — t'), 
the intensity predicted by the transport simulation 
for an instantaneous injection, with each Ii(t'), which 
yields the predicted response, Ri(t), due to each trian- 
gular injection (Figure 2b). Now we want to consider 
the response, R(t), to a general linear combination 
of the triangular injections. Because the transport 
equation is linear in F, the response to a linear com- 
bination of injections, I(t') — J2i a ih{t'), is the linear 
combination of responses, R(t) = ^aji?j(t). There- 
fore we can use linear least-squares fitting to find the 
linear combination that minimizes the x 2 between 
R(t) and the observed data (Figure 2d). Because each 
Ii(f) has a peak value of one, the coefficients, a,, are 
the values of the injection function at each joint t.- L 
(Figure 2c); the least squares fit also directly yields 
the statistical errors of these values. 

Initially, we set to to the peak time of the Ha flare, 
set t\ and £2 so that t^ — t\ and t\ — to were equal to 



the width of the time bins of the data, and set fur- 
ther ti so that each interval was twice the preceding 
interval. The joints, ti, were then adjusted according 
to an objective procedure. While the procedure could 
have been fully automated, we considered it prudent 
to manually execute fits for each set of joints, exam- 
ining each fit by eye; we never found it necessary to 
contradict the decisions mandated by this procedure. 
Each fit runs in the blink of an eye on a Pentium 
processor and the fitting procedure is completed in 
minutes. 

The second deconvolution technique assumes an in- 
jection function of the form 

I(t) = [C/(t - t Q )} exp[-A/(t - t ) - (t - t )/B], 

where A, B, C, and to are free parameters. This form 
was originally proposed by Reid [1964] as the solution 
of a coronal diffusion equation over the solar surface. 
However, we adopt this so-called Reid profile only as a 
convenient and widely understood parameterization, 
and we stress that we do not assume the existence of 
coronal diffusion. 

We perform a nonlinear least squares fit in which 
for each set of parameter values, A, B, and to, we 
numerically evaluate R(t) from (|J) using the results 
of the transport simulation, and calculate \ 2 for the 
fit to the data. (Given the other parameters, the op- 
timal value of C can readily be determined.) The 
parameters are optimized by the conjugate direction 
method [Press et al., 1988]. In practice, it was nec- 
essary to limit the variation of A, B, and to to phys- 
ically reasonable values. Unfortunately, for nonlinear 
optimization, one cannot be certain that the global 
minimum has been found, and we had to restart each 
fit several times for different initial parameter values. 
Fits were performed for various values of A so as to op- 
timize x 2 - Each fit required ~ 2 hours on a Pentium 
processor, making the entire procedure much slower 
than the piecewise linear technique. Although both 
deconvolution techniques yielded good fits and consis- 
tent results, we prefer the piecewise linear deconvolu- 
tion technique, which is faster, does not require any 
subjective evaluation by the researcher, and permits 
a more flexible shape for the injection function. 

3. Observations and Results 

We have examined data on protons stopping in the 
ISEE 3/MEH detector. During the times considered 
here, ISEE 3 was located near the inner Sun-Earth 
Lagrangian point. Two types of data have been used. 
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Priority rate (PR) data consist of raw counts in eight 
directional sectors over 96 s intervals that satisfy the 
priority 1 logic (Dl x D3 x D2 + D13A x^HxD5xD6 
[Kroeger, 1986]). It expected that during an intense 
solar event this rate is dominated by protons (other 
ions and electrons have priority 2), and detector simu- 
lations indicate that the stopping energies are roughly 
27-147 MeV. Large PR counting rates may have a 
very small statistical error, yet there are still fluc- 
tuations that are probably related to magnetic field 
irregularities. In that case, the statistical errors do 
not reflect the true level of effectively random rate 
fluctuations, leading to large \ 2 values even for a rea- 
sonable fit. Pulse height (PH) data enable one to 
reject background events more cleanly and also to de- 
termine the energy (within ~ 1 MeV) and time of 
arrival of each particle. While such data are clearly 
more desirable, the transmission of PH data to Earth 
was constrained by telemetry limits and depended on 
the priority logic. One must correct for this, and dur- 
ing times with a large electron flux, the PH data may 
yield poor statistics for the lower-priority protons. 

We chose to test the deconvolution techniques us- 
ing data for two solar events, one impulsive (X ray 
duration < 1 hour) and one gradual (X ray dura- 
tion > 1 hour). Thus the objective is to examine 
events for which the particles are evidently of solar 
origin, and are sufficiently intense to permit a detailed 
analysis and for which there is an evident anisotropy, 
which is an aspect of the fitting we would like to test. 
Other considerations are the absence of data gaps at 
the very start of the event and a small coronal dis- 
tance between flare site and the footpoint of the (av- 
erage) magnetic field connected to the spacecraft. We 
chose to examine the gradual event of July 20, 1981, 
and the impulsive event of January 2, 1982. Starting 
with such "well-connected" events simplifies the in- 
terpretation of the results; in future applications one 
could use these deconvolution techniques for poorly 
connected events in order to study lateral transport 
mechanisms. 

Figure 1 shows PR proton data for the gradual 
event of July 20, 1981. This event was associated 
with an Ha flare at 25°S, 75°W, peaking at 1322-1336 
UT (based on two observatories, Solar-Geophysical 
Data). The X ray decay time was 67 min [Oliver 
et al., 1989], indicating a gradual flare. The solar 
wind speed was approximately 375 km s^ 1 (S. Bame, 
private communication via ISEE 3 data pool, 1981). 
The magnetic field varied rapidly in both magnitude 
and direction during the time of interest (E. Smith, 



private communication via ISEE 3 data pool, 1981). 

Although the intensity was rather smooth as a 
function of time, the anisotropy times intensity sud- 
denly fell by a factor of about 4 at 1600 UT and re- 
covered at 1700 UT. Such sudden disappearances of 
anisotropy have been observed for other flare events 
by Evenson et al. [1982]. As discussed earlier, the 
anisotropy vector closely follows the magnetic field 
direction, so it would be expected to be especially 
sensitive to erratic magnetic field fluctuations. In 
fact, the field magnitude dropped particularly sharply 
during 1600-1700 UT, corresponding to a very short 
focusing length, L = -B/(dB/dz) w -0.04 AU; 
this strong reverse focusing apparently negated most 
of the outward-going anisotropy during that time. 
It is impossible for a transport model based on an 
Archimedean spiral field to predict such drastic fluc- 
tuations in the anisotropy. Therefore we conclude 
that when the magnetic field is erratic, as for this 
flare, our transport simulations are only appropriate 
for fitting the intensity as a function of time and par- 
ticle energy. 

For the proton event of July 20, 1981, there were 
sufficient statistics to use PH data and to sepa- 
rately examine different proton energy ranges. Fits 
were performed using both a piecewise linear injec- 
tion profile and a Reid profile, assuming a position- 
independent radial mean free path, X r , as recom- 
mended by Palmer [1982]. The resulting fit parame- 
ters are given in Tables 1 and 2, and piecewise linear 
injection functions are shown in Figure 3. In each 
case, both deconvolution techniques yielded the same 
best fit value of A r , except for the highest energy bin 
where they differed slightly. There is a hint of an 
increasing trend of A r with energy, which is consis- 
tent with previous results and theoretical calculations 
[e.g., Droge et al, 1997; Schmidt and Droge, 1997] 
The injection functions can be compared in terms of 
the full width at half maximum (FWHM). While the 
Reid profile always has a higher FWHM, the energy 
dependence is similar for both methods. 

Both deconvolution techniques indicate a much 
narrower injection profile for the highest energy bin 
(75-147 MeV). As discussed earlier, there is evidence 
that for gradual events the acceleration of ions takes 
place at a CME-driven shock as it propagates out- 
ward through the corona. Therefore the injection as 
a function of time can also be interpreted as injec- 
tion as a function of distance [Kahler et al, 1990], 
though the CME speed is not known for this partic- 
ular event. Thus these results suggest that for this 
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event the CME/shock system lost the ability to ac- 
celerate particles to ~100 MeV after traveling beyond 
a certain distance from the Sun. 

By integrating the injection function over time 
for each energy bin, we can estimate the spectrum 
of emitted particles. We compare our techniques 
with the commonly used time of maximum (TOM) 
method, in which one uses the observed peak inten- 
sity for each energy interval to estimate the relative 
spectrum. For the July 20, 1981, event, we find that 
the piecewise linear and Reid injection profiles yield 
very similar absolute spectra, except at the highest 
energies (Table 3). The relative spectra of both meth- 
ods are also similar to the TOM spectrum, except 
at the high energy bin. The TOM technique implic- 
itly assumes a similar injection profile and scattering 
mean-free path for each energy, so one might well ex- 
pect some deviation for the high energy bin, where 
both deconvolution techniques yield a much shorter 
duration of injection. Also, there were large statis- 
tical errors for this energy bin. From our results we 
conclude that the TOM relative spectrum can be ac- 
curate for a short energy span (here a factor of 3), 
but for wider energy ranges it is worthwhile to deter- 
mine the spectrum more accurately by deriving the 
energy-dependent injection profiles. Note that the 
TOM method does not yield an absolute spectrum, 
and the consistency of our two methods for the ab- 
solute spectrum improves our confidence in both of 
them. 

We have also analyzed PR data from the impul- 
sive flare event of January 2, 1982 (Figure 4). This 
event was associated with an Ha flare at 18°S, 88° W, 
peaking at 0620-0621 UT (Solar-Geophysical Data). 
The X ray decay time was 16 min [Oliver et at, 1989], 
and the solar wind speed was approximately 350 km 
s _1 (S. Bame, private communication via ISEE 3 
data pool, 1982). The magnetic field was stable in 
magnitude and direction during the time of interest 
(E. Smith, private communication via ISEE 3 data 
pool, 1982). Because of this, we were able to analyze 
the anisotropy in addition to the intensity for this 
event. Unfortunately, because of the limited statis- 
tics and strong intensity of electrons (which had a 
higher priority for pulse-height telemetry), we were 
unable to examine PH data in detail to determine 
energy-dependent injection functions. 

Simultaneously fitting the intensity and anisotropy 
times intensity provides a stringent test of the injec- 
tion and transport models. Our two techniques for 
deconvolving the effects of interplanetary transport 



yielded similar FWHM durations of injection (17 and 
24 min, respectively; see Tables 1 and 2) and yielded 
reasonable fits to the intensity and anisotropy times 
intensity. The duration of injection over this broad 
energy range (27-147 MeV) was similar to the decay 
time of X ray emission (16 min) and was markedly 
shorter than the corresponding duration for the grad- 
ual event of July 20, 1981, though that event had a 
similarly short injection duration for the highest en- 
ergies. 

We close this section with examples of how less ac- 
curate transport assumptions can artificially broaden 
the derived injection profile. For an analogy, consider 
the deconvolution of a telescope image to account for 
the point spread function. An inaccurate estimate of 
the point spread function yields a deconvolved image 
that is still artificially broad compared to the true im- 
age size. When an improved point spread function is 
used, one obtains a sharper image. 

Figure 5 shows examples of a similar effect for our 
deconvolution problem, in which we compare fits to 
PR data from both events, assuming that either A r 
or Am is constant in position. Previous authors [e.g., 
Palmer, 1982] have concluded that a constant A r is 
a better and reasonable assumption for transport in 
the inner heliosphere. We find that the fits assum- 
ing a constant A|| (Figures 5c and 5d) yield much 
broader injection functions than those based on the 
presumably more accurate assumption of a constant 
A r (Figures 5a and 5b). These examples stress the im- 
portance of using an accurate transport model when 
determining the injection function near the Sun. 

4. Discussion 

The previous study with goals most similar to ours 
was that of Ma Sung and Earl [1978], which em- 
ployed approximate analytic solutions to a focused 
transport equation. One of their assumptions was 
a position-independent ratio, Am/L, where L is the 
focusing length. This assumption was necessary for 
their analytic approximation but is less accurate than 
the assumption of a constant A r and a function L(z) 
derived for an Archimedean field. Based on our re- 
sults, one would expect injection functions that are 
smeared out when compared with those for a con- 
stant A r (Figure 5). Those authors frequently found 
particle release times of the order of hours, and it is 
possible that their injection profiles were artificially 
broad due to that assumption. The deconvolution 
techniques presented here, along with more realistic 
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transport assumptions and modern numerical simu- 
lations, could be profitably applied to survey the in- 
jection profiles of a variety of events, as Ma Sung and 
Earl [1978] did. Shortly, we will consider how such 
results should be interpreted in light of our modern 
understanding of the origin of solar energetic parti- 
cles. 

A valid question is how one knows that injection 
occurs near the Sun, when the interplanetary shocks 
that frequently accompany gradual events are known 
to be capable of accelerating particles [Gosling et al., 
1981]. For the gradual event considered here, several 
days after the flare, the log(proton intensity) versus 
time showed a bump (July 23) and then a double- 
peaked increase (July 24-25) to fluxes <l/30 of the 
initial peak. There was intense geomagnetic activity 
at these times (Solar- Geophysical Data). Presumably 
these features are associated with the passage of a 
CME/shock system, which affected the particle prop- 
agation or accelerated particles to higher energies. In 
this work, we focus on the proton observations on 
July 20-21, 1981, well before such features arose. The 
rapid rise and exponential decline are consistent with 
injection shortly after the flare occurrence convolved 
with the effects of interplanetary transport. Since the 
CME could not have traveled far from the Sun during 
that time, we concluded that the emission was near 
the Sun and neglected any source motion. 

We note further that all measurements of ionic 
charge states above 3 McV/nucleon indicate that so- 
lar energetic ions are accelerated out of coronal ma- 
terial [Boberg et al., 1996]. This implies that the 
CME/shock system first accelerates ions out of coro- 
nal material, presumably while it is still near the Sun, 
and can also further accelerate ions as it propagates 
outward. The relative importance of acceleration at 
different distances from the Sun undoubtedly varies 
from event to event and varies with energy. 

Our results indicate that for the gradual event of 
July 20, 1981, the duration of injection was much 
shorter for higher energies (Figure 3), which implies 
that the CME/shock system lost the ability to ac- 
celerate a significant flux of ~100 MeV protons after 
traveling a certain distance from the Sun. We also 
note that the shock-associated particle increases on 
July 23, and July 24-25 had a significantly steeper 
spectrum than the main peak on July 20 and were not 
seen at all in the high energy bin. Therefore one might 
interpret the short duration of injection at higher en- 
ergy as indicating a transition from injection with a 
harder spectrum (when the CME shock is close to the 



Sun) to injection with a softer spectrum (as for the 
interplanetary shock). While it is imprudent to draw 
general conclusions from one flare/CME event, con- 
firmation of this result for other events would provide 
important information on the acceleration mechanism 
for escaping energetic ions from such events. 

It is also interesting to compare the best fit val- 
ues of A for fits assuming a constant Am or a con- 
stant A r , which are largely determined by the ob- 
served intensity decay at long times. We obtain ratios 
of A||/A r = 5.8 and 4.5 (Table 1). At any given point, 
X\\/X r = sec 2 ip, so the observed ratios are charac- 
teristic of r — 1.9 and 1.5 AU. This indicates that 
the intensity versus time of solar particles at 1 AU 
is strongly influenced by the transport conditions at 
greater distances. 

Since we have only treated the interplanetary trans- 
port parallel to the magnetic field, we have derived 
injection profiles under the assumption that parti- 
cles are strictly confined to a narrow flux tube con- 
necting the observer to the source near the Sun. 
However, perpendicular diffusion implies that the ob- 
served particles could have originated from a distribu- 
tion of longitudes and latitudes. Summarizing a vari- 
ety of results, Palmer [1982] recommends using k± = 
(v/c) x 10 21 cm 2 s _1 at 1 AU (well away from sector 
boundaries, corotating interaction regions, etc.). For 
our typical proton velocity of c/3, this implies an an- 
gular spread of ~7° after 3 hour or ~20° after 1 day. 
In comparison, gradient and curvature drifts imply an 
angular motion ^0.1 AU day -1 near the Earth or an 
integrated angular motion ~0.1° for the first particles 
that arrive and are hence negligible during the time 
of interest, as is the rotation of the Sun. 

For gradual events, there is evidence that coronal 
mass ejection shocks can accelerate and release parti- 
cles over a large fraction of the solar surface [Mason 
et al, 1984; Kahler, 1992], so our results are really 
telling us the injection profile averaged over ~ 7° of 
solar latitude and longitude from the footpoint mag- 
netically connected to the observer (which is itself un- 
certain, except when there are observations of mov- 
ing interplanetary type III radio bursts that track 
the mean motion of electrons along the interplane- 
tary magnetic field [Reames and Stone, 1986; Reiner 
et al., 1995]). The accuracy of the measured dura- 
tion of injection as a function of energy should be 
unaffected insofar as perpendicular diffusion does not 
significantly increase the distance traveled before par- 
ticles actually reach the observer, i.e., as long as it is 
much weaker than parallel diffusion. 
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For impulsive events, the interpretation is po- 
tentially more complex: Since particles are acceler- 
ated at the flare site, their lateral spread could be 
due to perpendicular diffusion in the interplanetary 
magnetic field (including the "random walk" of the 
magnetic field lines themselves [e.g., Jokipii, 1966; 
Matthaeus et at, 1995]) or to lateral transport within 
the corona. For impulsive flares that are magnetically 
well-connected to the observer, such as the event of 
January 2, 1982, which we have considered here, we 
can again argue that our determination of the dura- 
tion of injection is accurate as long as perpendicular 
diffusion does not significantly increase the propaga- 
tion time. However, the absolute normalization of the 
observed flux depends strongly on the lateral extent 
of the particle distribution at the observer's radius. In 
this case, we are determining the number of particles 
injected per solid angle at the observer's radius. 

5. Conclusions 

We have developed two techniques for deconvolving 
the effects of interplanetary transport using numeri- 
cal solutions of the transport equation: (1) assum- 
ing a piecewise linear injection profile and (2) using 
a Reid injection profile. The deconvolution can yield 
the interplanetary scattering mean free path, the in- 
jection profile as a function of energy, and the spec- 
trum. The two techniques yield consistent results for 
the gradual flare/CME event of July 20, 1981, and 
the impulsive flare event of January 2, 1982, giving 
us confidence in both techniques. The piecewise lin- 
ear profile is preferred because the fitting procedure 
is faster, objective (not relying on user evaluation of 
fits), and permits a more flexible profile shape. It is 
important to examine the observed magnetic field; if 
this is erratic, it may not be possible to fit the ani- 
sotropy data with simple transport assumptions. For 
the July 20, 1981, event, a simple time-of-maximum 
estimate of the spectrum agreed with our results over 
a factor of 3 in energy but deviated at higher ener- 
gies. The FWHM of injection was much shorter for 
higher energies (~100 MeV), indicating that for this 
event, the CME and associated shock lost their effi- 
ciency of accelerating such particles after traveling a 
certain distance from the Sun. 
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Figure 1. Intensity and anisotropy times intensity as a function of time for 27-147 MeV protons following the 
July 20, 1981, gradual solar flare event. 

Figure 2. Illustration of the deconvolution technique for a piecewise linear injection function near the Sun. The 
transport equation is solved for an instantaneous injection of particles. The resulting Green's function is convolved 
with (a) triangular injection profiles to (b) yield response functions. Linear, least squares fitting yields (d) the 
linear combination of response functions that best fits the data and (c) the corresponding best-fit piecewise linear 
injection profile. 

Figure 3. Fits to the (a) observed proton intensity versus time in four energy bins for (b) optimal piecewise 
linear injection profiles. Note the expanded timescale in the right panels. 

Figure 4. Fits to the (a) observed 27-147 MeV proton intensity versus time and anisotropy times intensity versus 
time for the (b) optimal piecewise linear injection profile. Note the expanded timescale in Figure 4b. 

Figure 5. Comparison of best-fit piecewise linear injection profiles for transport simulations assuming either 
(a-b) A r or (c-d) A|| to be independent of position for 27-147 MeV protons on July 20, 1981, for Figures 5a and 5c 
and January 2, 1982, for Figures 5b and 5d. The assumption of a constant A r is expected to be more accurate, and 
it yields a sharper injection function for both events. 
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Table 1. Fit Parameters for a Piecewise Linear Injection Profile 



Date 


Data 


Energy, 


A r or Am 


A, 




X 2 /di. 


FWHM, 


Figure 






MeV 


constant? 


AU 






min 




July 20, 1981 


PH 


27-39 


\ r 


0.08 


{7.2,0.5,0.7} 


1.92 


33 


3 




PH 


39-55 


X r 


0.10 


{1.6,1.9} 


0.96 


65 


3 




PH 


55-75 


\ r 


0.12 


{0.57,0.93} 


1.11 


52 


3 




PH 


75-147 


\ r 


0.10 


{0.58,0.04} 


1.24 


12 


3 




PR 
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\ r 
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2, 5 
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89.68 
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{12.0,2.7}xl0" 3 
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4, 5 




PR 


27-147 
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2.35 
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a Units of 10"° sr" 1 s" 1 MeV' 1 . 



Table 2. Fit Parameters for a Reid Injection Profile 



Date 


Data 


Energy, 
MeV 


A r , 
AU 


A, 
min 


B, 
min 


C a 


to, 
UT 


X 2 /d.f. 


FWHM, 
min 
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24 



a Units of 10^ sr" 1 s" 1 MeV" 1 . 



Table 3. Spectra of Protons Injected on July 20, 1981 



Technique Energy Range, MeV 

27-39 39-55 55-75 75-147 



Absolute Spectra, 10 26 sr" 1 s" 1 MeV" 1 

Piecewise linear 19,000 7120 2770 465 
Reid profile 19,900 7120 2810 581 

Relative Spectra 



Piecewise linear 1 0.374 0.145 0.0244 

Reid profile 1 0.358 0.141 0.0292 

Time of maximum 1 0.343 0.136 0.0185 
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